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Abstract. Baxter's TQ-equation is solved for the six-vertex model using the 
representation theory of quantum groups at roots of unity. A novel simplified 
construction of the Q-operator is given depending on a new free parameter. 
Specializing this general construction to even roots of unity and lattices with an 
odd number of columns two linearly independent operator solutions of Baxter's 
TQ equation are obtained as special limits. 



1. Introduction 

Integrable lattice models of statistical mechanics form an important interface 
between physics and mathematics and have given rise to numerous discoveries relat- 
ing both subjects, quantum groups are a celebrated example |Dr86ll3i85| . Recent 
discoveries of infinite-dimensional non-abelian symmetries at roots of unity in the 
case of the six-vertex model |DFM04] and new developments in the eight-vertex 
model FM03, FM04al IFM04b] have opened up a number of questions of math- 
ematical and physical interest. The degeneracies in the spectrum of the transfer 
matrix at these particular points and their relation with Bethe's ansatz |Be31| and 
Baxter's TQ-equation |Ba72l lBa73l IB"a82j have been subject of a series of papers 
FMOla, FMOTBl lB"a02l IBa04] . Moreover, the case of primitive roots of unity of 
order three has found particular interest because of its connection with combinato- 
rial aspects such as the enumeration of alternating-sign matrices or plane partitions, 
see e.g. |Ku96llK50ll BdGNOlj. Here Baxter's TQ-equation has been explicitly 
solved for the groundstate eigenvalue of the six-vertex [FSZOOL IStOl) and more 
recently also for the eight- vertex model |FM04bl IBM04| . 

In an independent development the techniques of the quantum inverse scat- 
tering method FST79 have been applied to conformal held theory and Baxter's 
TQ equation has been discussed in the context of the Liouville model BLZ99, 
FKV01 . This led to a representation theoretic approach of constructing Baxter's 
Q-operator and in subsequent papers the method has been generalized and applied 
to the six-vertex model away from a root of unity |AF97L IRW02L lKo04b] . 

The representation theoretic construction of Q-operators for the six-vertex 
model at roots of unity and their relation to the aforementioned symmetries has 
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therein that the Q-operators constructed in |Ko03| provide an efficient tool to 
analyze the degeneracies in the spectrum of the transfer matrix, reveal representa- 
tion theoretic information on the affine symmetry algebra and allow to determine 
whether a second linear independent solution to Baxter's TQ equation exists. It is 
the last aspect which we will further highlight in this article. 

The possible existence of two linearly independent solutions to the six-vertex 
model has been discussed in |PS99j excluding the root of unity case. See also 

KLWZ97 for the discuss ion of the eight- vertex case. While the aforementioned 
papers addressed the question on the level of eigenvalues, it has been shown in 
|Ko04c| that these solutions arise in the spectrum of the Q-operators explicitly 
constructed in |Ko03| at primitive roots of unity. The analogous discussion for 
the eight-vertex model based on Baxter's Q-operator constructed in his 1972 work 

Ba72 and numerical computations has been subject of the paper FM04B]. There 
is, however, a significant difference as Baxter's construction procedure for the Q- 
operator does not rely on representation theory or contains free parameters. It 
is the representation theory of quantum groups at roots of unity used in |Ko03] 
which provides the key to the derivation of the TQ equation and a second equally 
important functional equation, see equation 1|5.1[) in the text. 

In comparison with the discussion contained in |Ko04c| the new aspect pre- 
sented here is the explicit construction of Q-operators which contain one more free 
parameter. This puts one into the position to take two special limits which yield the 
two linearly independent operator-solutions of Baxter's TQ equation; see Definition 
3.2 and Definition 3.5 in the text. Based on these results we will show for even roots 
of unity and an odd number of lattice columns that the proof of existence of two 
linear independent operator-solutions to the TQ equation and the completeness of 
the Bethe ansatz above and below the equator can be reduced to a well-defined 
mathematical problem: the construction of an intertwiner for two quantum group 
representations at roots of unity. In the appendix we then explicitly construct this 
intertwiner for the cases where the order of the root of unity is four and six. 



Our definition of the six-vertex model and the discussion of Baxter's TQ- 
equation will be based on the Drinfel'd-Jimbo quantum group associated with the 
affinc algebra s?2- We briefly recall some of its defining relations in terms of the 
Chevalley-Serre basis; see e.g. CP94 for further details. 

Proposition 1. There is a quasi-triangular Hopf algebra U q (sl2) which is gen- 
erated by elements {e^, fi, q ±hi }j=o,i obeying the relations 



2. Preliminaries 
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^See also IBS9 for a different construction and the related comments in the introduction 
of |KoOS| . 



SOLVING BAXTER'S TQ-EQUATION 



3 



with Aij = (— iy +J 2 being the Cartan matrix. In addition, fori ^ j the q-deformed 
Chev alley- S err e relations hold, 

(2.4) x^Xj [SjgX^x^x^ -\- \3^qXiXjX^ XjX.^ — 0, x — e, f . 
The coproduct ofU q (sl<2) is given by 

(2.5) A(e;) = l®ei + q hi ®ei, A(/0 = ft ® q- h > + 1 A(q hi ) = q hi ®q hi . 

The opposite coproduct A op is obtained by permuting the two factors. There exists 
an universal R-matrix intertwining these two coproduct structures 

(2.6) RA(je) = A op (x) R, x£U q (d 2 ), R £ U q (b + ) ® U q {b-) . 
Here U q (b±) denote the upper and lower Borel subalgebra, respectively. 

Before we can define the six-vertex model in terms of representations of the 
quantum group U q {sl2) we need one more ingredient, Jimbo's evaluation homo- 
morphism. 

Proposition 2. Let z £ C be nonzero then the mapping 

(2.7) e - z f, f Q - z- l e, q h ° -> if"*, e x -> e, /i - /, - q h . 

defines an algebra homomorphism ev z : Uq^slz) — > U q (sl2)- Here U q (sl2) is isomor- 
phic to the Hopf algebra generated by either {ei, fx, q hl } or {eo, /o, q h °}. 

Denote by 7ri"' ) = o eu z : Uqislz) — > EndC ,l+1 the spin n/2 evaluation 
representation of the quantum group defined by the relations 

7r (n) (e)|m) = [n-m+l] q \m-l), 

(2.8) n-W (/) |m) = [m + 1], |m + 1) , tt^ (g' 1 ) |m> = q n ~ 2m \m) . 

Here {|to)}J^ =0 denotes the standard orthonormal basis in C™ +1 . For each integer 
n > consider the intertwiner of the tensor product iri^ (g) tt^, i.e. 

(2.9) i?(" +1 )(z) = (ttW (8 jrf 5 ) R € End (C" +1 ® C 2 ) , 

whose matrix elements w.r.t. the second factor are calculated to be (up to an overall 
normalizing factor) 

(0|i?(" +1 ^) |0) = zq^ n \q h ' 2 )-^ n \q- h ' 2 ), 

(01^"+%) |1) = «?(«-g-V n V / V n) (/), 
(l|i^+%)|0) = (g- g - 1 ) 7 rW(e)7rW(g- h / 2 ), 

(2.10) (l|i? (n+1) (z)|l) = zq^ n \q- h / 2 ) ~^ n \q h ' 2 ) . 
We now introduce the six- vertex fusion hierarchy K R87) . 

Definition 2.1. Define the six- vertex fusion matrix T^ : (C 2 )' 8M — » (C 2 )® 11 
of degree n £ N by setting for some fixed integer M > 1 

(2.11) T^\z) = Tr R^(m n )---R&\z<l n ) ■ 
The two special elements 

(2.12) T^\zq- 2 ) = T{z) and T«(z) = (zq 2 - 1) M id , 

are called the six- vertex row-to- row transfer matrix T and the quantum determinant 
KRS81 TW, respectively. 
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PROPOSITION 3 ( |KR87| ). The fusion matrices satisfy the functional equation 

(2.13) T^\z)T^ 2 \zq- 2 ) = {zq 2 - l) M (zq~ 2 ) + (z - 1) M T^~ l \zq 2 ) 

and hence can be successively generated from the transfer matrix and quantum de- 
terminant. 

Proof. The assertion follows from the exact sequence |CP94 , 

(2.14) -> tt^- 13 ^ ttW ® 7T« ^ tt^, + 1) -> 0, «; = w'q- 1 = w" q = zq n+1 

See also Ko04b for details such as the explicit form of the inclusion i and projec- 
tion map p needed to derive l|2.13|l . □ 

The statistical mechanics model is now defined in terms of the partition function 
Zq v = TrT M which is the physical quantity of interest. Here M is the number of 
lattice columns and M' the number of lattice rows. In this definition we have 
imposed periodic boundary conditions on the lattice. This model is said to be 
integrable as the transfer matrix T has infinitely many conserved quantities. 

Corollary 1. From the definition $2.11\) and the existence of the universal 
R-matrix it immediately follows that 

(2.15) [T ( - m \z),T ( - n \w)] =0, Vz, meC, m,neN. 

Proof. The assertion is a direct consequence of the Yang-Baxter equation, 

R.I2R13R23 = R.23R-13R-12- n 

For later purposes we note further symmetries of the six- vertex fusion matrices. 
Their proof follows by direct calculation. 

Lemma 2.2. Let o~ x ,o~ y ,o~ z denote the Pauli matrices. Define the following 
matrices acting on (C 2 )® M ; 

M M M 

(2.16) s* = - « = n <' 6 = n < ■ 

m— 1 m— 1 ra— 1 

Then for any integer n> 1 the fusion matrix T^ commutes with S z ,${ and 6. 

Throughout this article we will implicitly assume that we always work in a basis 
of eigenvectors of the total-spin operator S z , i.e. we will treat S z as a diagonal 
matrix with half-integer entries in the interval [—M/2,M/2]. 

Lemma 2.3. The transfer matrix satisfies the identity 

(2.17) T(z,q) = T(zq 2 ,q- 1 ) t 
and hence is normal provided \q\ — 1; see $2.15\) . 

3. The Auxiliary Matrix 

Henceforth we will specialize to the case when q is a primitive root of unity of 
order N > 3. The six-vertex fusion hierarchy stays well-defined in this limit. We 
set N' = N if the order is odd and N 1 — N/2 if it is even. 
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Proposition 4. For any r, s, z € C the following defines an N' -dimensional 
representation n + = 7r + (z; r, s) of the upper Borel subalgebra U q (b+): 

n+(q hl )\n) = TT+(q-' 10 ) \n) = rq- 2n \n) , 
7r+(e )H = z\n + l), 7r+(e )|iV'-l)=0, 

„ I 1 _ „2n _ S(7 -2n 

(3.1) 7r + ( ei )|n) = (g-g- 1 ) 3 ' n ~ 1) ' ^( ei )l°) =0 ' 

-ffere {l")}^^ 1 denotes the standard orthonormal basis in <C N . Let lo denote the 
algebra automorphism {ei, eo, q , q h °} — ► {eo, ei, g' i0 , q hl } and set 

(3.2) 7T~ := 7T + O . 

A short calculation shows that the defining relations of the quantum group are 
satisfied in the representations tt^. In order to make contact with previous results 
in the literature we have the following 

Remark 3.1. Note that the representation ir~ is a particular root-of-unity 
restriction of the representation in RW02 . For the choice of parameters r — 
H~ 1 q~ 1 and s = /i~ 2 the representation l|3.1|l coincides with the representation 
used in |Ko03l IKo04al IKo04bl IRo04c] . 

Let the matrix L = a® a + a~ + [3 <E> a + + 7 ® o~ + S (g> (T~a + be the intertwiner 
of the tensor product ir + ® tt^2i of evaluation representations, explicitly 

L (tt + ® ?rf 3 )A(a;) = ((vr+ ® 7rf 3 )A° P (z)) L, Va; € U q (b+) . 

From this relation the matrix elements up to an overall normalization factor are 
computed to 

a = zs/rn + (q^) - n + (q~^), (3 = (q - q^ 1 )^ (e )TT + (q~ ^ ) , 

(3.3) 7 = (q-q-^n+ie^n+iq-^), 6 = z rq 2 tt+ (<T* ) - tt+ (q^ ) . 

Notice that (up to a simple gauge transformation) the intertwiner for the rep- 
resentation 7r~ is obtained via the permutation {a, f3, 7, 8} — > {8, 7, f3, a} which 
corresponds to spin reversal, i.e. L — ► (1 ® a x )L(l <£> a 1 ). 

Definition 3.2. Define the auxiliary matrix in terms of the intertwiner L as 
the trace of the following operator product, 

(3.4) Q(z:r,s) = Tr L 0M (z;r,s) ■ ■ ■ L 01 (z;r,s) . 

This matrix commutes by construction with the fusion matrices as the inter- 
twiner L must satisfy the Yang-Baxter equation DJMM90 

L 12 (z/w)L 13 (z)R%{w) = Rg\w)L 13 (z)L 12 (z/w) . 
Furthermore, Q preserves two of the symmetries l|2.15ll Ko03, Ko04a j , i.e. 

(3.5) [Q(z; r, s), T^ n \w)} = [Q(z; r, s),S z ] = [Q(z; r, s),&) = . 

Spin-reversal symmetry on the other hand is broken |Ko03L IKo04aj . To see this 
we first note that dependence of the auxiliary matrix Q on the parameter r can be 
easily extracted. 

Lemma 3.3. The auxiliary matrix can be simplified as follows 

(3.6) Q{z- r, s) = r~ s 'Q(z; 1, s) = r- 3 *Q(z; s) . 
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Proof. Because the auxiliary matrix preserves the total spin a general matrix 
element {<r\, ...,<tm\Q\qi, Qm) must contain as many /J-matrices as it does 7- 
matrices whose r-dependence cancels against each other. Each occurrence of an a- 
matrix contributes a factor r" 1 / 2 while each <5-matrix gives rise to a factor r 1 / 2 . □ 

As an immediate consequence of the decomposition of the auxiliary matrix 
w.r.t. to the parameters r, s we have the following 

COROLLARY 2. The auxiliary matrix constructed in this article can be identified 
with the auxiliary matrix in |Ko03L IKo04a( IKo04c| via the relation 

(3.7) Q^(z) = f i s *q s *Q(z;^ 2 ) ■ 
We now easily find the identities 

Lemma 3.4. Under spin reversal the auxiliary matrix transforms as 

mQ(z,q; S )m = (-z^qM+^s^Qiz-^s-^q^y 

(3.8) = Q{zq 2 s,q- 1 ;s- l ) t = q 2S 'Q{zs,q ]S - 1 ) . 

From the last equation the conjugate transpose of the auxiliary matrix is deduced to 
be 

(3.9) Q(z, q; s)* = Q(z, g" 1 ; sf = q 2S *Q(zq~\ q; s) . 

Proof. The transformation properties l|3.8|) are a direct consequence of the 
identities in |Ko03L IKo04al IKo04cj and ftTTjL □ 

As we want to solve the eigenvalue problem of the six- vertex transfer matrix in 
terms of the auxiliary matrix we need to show that Q is diagonalizable. Further, 
we wish to establish the analytic properties of the eigenvalues. 

Claim 1. Let z,w,s,t 6 C be arbitrary independent complex numbers. Then 
the following commutation relation holds for all primitive roots of unity 

(3.10) [Q(z;s),Q(w;t)] =0 . 

Proof for N = 3,4, 6. A sufficient (not necessary) condition for the above 
commutation relation to hold is to show that the intertwiner, say S 1 , for the tensor 
product 7r + (z; 1, s) ® tt + (w; l,i) exists. If it does, S is bound to satisfy the Yang- 
Baxter equation [DJMM90 Si2Li 3 (z; s)L 2 3(w; t) = L 2 3(w; i)Li 3 (z; s)S\2 from 
which the assertion immediately follows. Unlike the case when q is not a root of 
unity the existence of such an intertwiner cannot be deduced from the existence 
of the universal R-matrix; see CPQ4] for details and |Ko03| for the connection 
with the present discussion. It appears that at the moment the only way known 
to prove existence of S is through explicit construction. This has been carried out 
for the case N = 3 in Ko()4a . Here we state the intertwiners for N — 4, 6 in the 
appendix of this article. Numerical checks carried out for N = 5, 7, 8 and M < 11 
confirm (|3.10(l also in these cases. □ 

Henceforth, we take (|3.10|l as a working hypothesis. The two immediate impli- 
cations are: 

1. According to (j3.9|l the auxiliary matrix is normal, [Q(z; s),Q(z; s)*} — 0, and 
hence diagonalizable. 
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2. The eigenvalues of Q(z; s) are polynomials in z whose degree is at most M. By 
construction we know that (|3.4|1 can be decomposed as follows, 

M 

Qi z ,q; s ) = ^2 Q m ( s '«) zI "' 

m=0 

and l|3.10[l then ensures that the coefficients {Q m (s,q)} can be simultaneously di- 
agonalized. In addition we infer from Ij3.8(l that the coefficients are related via 

(3.11) Q M -m(s,q) = (-l) M g M - 2S %^ + ^g m ( S - 1 ,g- 1 ) . 

Finally, we can make use of the fact that the intertwiner i|3.3H is a lower-triangular 
matrix at z = from which we deduce the following identity for the zeroth term of 
the polynomial, 

N'-l 

(3.12) Q(0; s) = Q m=0 (s, <?) = (-l) M Tr n+ q- h * S * = (-1) M £ q 2lS * . 

1=0 

Note that because of l|3.11[l this expression determines at the same time whether 
Q has the maximal degree M or not. 

3.1. The limit s — ► of the auxiliary matrix. For later purposes we end 
this section by introducing two special limits of the auxiliary matrix. 

Definition 3.5. Let 7r ± = 7r ± (z;l,s) be the representations Q3.1|) . (|3.2|l and 
L the respective intertwiners ()3.3|l . Then the limit s — > is well-defined and we 
set 

(3.13) Q ± (z) := lim Tr n ±L± M {z) ■ ■ ■ L±(z) . 

From the remarks after equation l|3.3|l it follows immediately that the two 
operators are related via spin-reversal UiQ~ = The motivation for the above 

definition will become apparent in Section 5. 

4. Baxter's TQ Equation 

We are now ready to relate our construction of the auxiliary matrix to the 
six-vertex fusion hierarchy. The link is provided by a functional equation similar 
to Baxter's famous TQ equation. In contrast to his approach we can now exploit 
the representation theory of quantum groups in its derivation. 

Proposition 5. The auxiliary matrix j3.4\ ) and the six-vertex transfer matrix 
\2.1ty) obey the following operator q-difference equation for any primitive root of 
unity q, 

(4.1) Q(z; s)T(z) = (z - l) M q S *Q{zq 2 - sq' 2 ) + (zq 2 - l) M q~ S * Q{zq- 2 ; so 2 ) . 
Here S z is the total- spin operator defined in \2. lb)) . 

Proof. The strategy of the proof follows closely the one in Ko03 by means 
of the identification l|3.7|l . For the derivation of l|4.1|) it is more convenient to re- 
introduce the parameter r of the representation (|3.1|l and then at the end use (13. Gl) 
for taking the limit r — ► 1. One verifies the following non-split exact sequence of 
representations of U q (b + ), 

— > n + (zq 2 ; rq^ 1 , sq~ 2 ) A 7r + (z; r, s) (g) tt^ A n + (zq~ 2 \ rq, sq 2 ) — > . 
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This determines the functional equation by the same arguments as presented in 

[Ksnu. □ 

Notice that the above functional relation (|4.1() differs from Baxter's original TQ 
equation as his matrix does not depend on extra parameters which shift. However, 
as a trivial consequence we now deduce that the two auxiliary matrices l|3.13|l obey 
a functional equation which is of Baxter's type. 

Corollary 3 (Baxter's TQ-equation). In the limit s — > we have 

(4.2) Q ± (z)T(z) = (z- l) M q ±SZ Q ± (zq 2 ) + (zq 2 - l) M q^ S ' Q ± (zq' 2 ) . 

Notice that the equation for Q~ follows from the one for Q + by spin-reversal. 

Since transfer and auxiliary matrix are both simultaneously diagonalizable we 
can interpret the functional equations (|4.1|l . (|4.2() on the level of eigenvalues. 

Corollary 4. Let Q(z; s) be an eigenvalue of the auxiliary matrix \'d-4\) which 
is not identical zero. (We denote eigenvalues and operators by the same symbol.) 
Then the corresponding eigenvalue of the fusion matrix of degree n is given by the 
formula 

(4.3) *■>(,) - q M^, QHz)QH ^ , ± ^ ± . 

Proof. From (|2.13(1 the spectrum of the fusion matrices H2.11|l is calculated 
via induction. □ 

5. Another Functional Equation 

To fully analyze the structure of the eigenvalues of the auxiliary matrix Baxter's 
TQ-equation is not sufficient. Instead we make use of another functional equation 
which by the same philosophy as before follows from the decomposition of a tensor 
product of representations. 

Proposition 6. The auxiliary matrix \3-4\) obeys the functional equation 
(5.1) 

Q(zq 2 /s- s)Q( Z ;t) = Q(zq 2 /s- stq' 2 ) [{zq 2 - \) M + ^'(M-S')^'-!) ^3) 
for arbitrary complex numbers s,t £ C. 

Proof. The proof follows the same strategy as the one for the functional 
equation (|4.1|l . We invoke the identity (|3.7() and refer the reader to |Ko04c| for 

the derivation of the exact sequence 



-> 7r+(z'; s'iq' 1 , s') -> 7r+(zq 2 /s; s^q' 1 ^) ® ir + (z; t^q^ 1 ^) -» 

7r+(z"; S "^- 1 , S ")®^ + 2 1 ) -0 

where with regard to equation (51) in Ko04c we have set 

s = ^r 2 , t = v- 2 , s'i = (/^g)~\ s"i = (^vq- N '+ 1 )- 1 , z' = z" = zq 2 / s . 

The remaining proof now follows the analogous steps as in |Ko04c| . see the appen- 
dix therein, employing at the end l|3.7|l to arrive at (|5.1I) . □ 
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Corollary 5. Let M g 2N + 1 and N e 2N. Then all the eigenvalues of the 
auxiliary matrix are of the form 

-i 2N'S Z n M—n , r 

(5.2) Q{z;s) = - ~_ g Y[(l-zx+) Ylil-zsxj), 



y 4=1 4=1 
where the parameters x^ do only depend on q. This implies the following decom- 
position of the auxiliary matrix, Q(z;s) = Q~(0) _1 Q + {z)Q~(zs). Moreover, we 
have the identity 

N ' „-US z ( Zf] 2i _ i\M 

(5.3) Q-( Z ) = J* »Q-{0mz) g ^ { J^ +izq l 2) • 

Proof. We state a slightly simpler version of the proof given in |Ko04c| which 
applies to the case of arbitrary positive integers N,M. As N is even we have 
q N = — 1 and since M is odd the eigenvalues of the total-spin are half-odd integers, 
i.e. S z = -M/2, -M/2+1, ...,M/2. Thus, we can conclude from fTTfy and (|3~TT|) 
that the eigenvalues of the auxiliary matrix are polynomials in the spectral variable 
z of degree M and are non-vanishing at the origin. Hence, they can be written in 
the form 

n M — n 2N'S Z 

Q{.z-s)=M\{{l-zx+) H(l-z yj (s)), N = ~_ g 2SZ . 

4=1 4=1 q 

where the normalization constant Af = Q + (Q) = q~ 2S Q~{0) is given by l|3.12|l . 
(13.81) and l|3.13|l . Here we have split the zeroes of the eigenvalues into two groups: 
the zeroes x^~ do not depend on the parameter s, while the zeroes yi = yi(s) do 
depend on it. We allow for the extreme cases n = and n = M. Employing the 
functional equation H5.lt we find that the ratio 



Q{zq 2 /s;s)Q(z;t) _ ^ ^ {I - zq 2 V] {s) / s){\ - zy 3 {t)) 

4=1 4= 



Q(zq 2 /s;stq 2 ) ' fj- 3 fj- I - zq 2 yj{stq 2 )/.s 



must be independent of the arbitrary complex parameters s, t which implies that 
the zeroes yi depend linearly on s, i.e. j/j(s) = x~ s for some x~ which only depends 
on q. Invoking now once more the transformations (|3.8|) to determine n this proves 
the first assertion. The second now trivially follows also from (|3.8|) and l|3.13|l . 
The third and last identity is immediate from the functional equation l|5.1|l and the 
expression (14.31) for the spectrum of the fusion matrices. □ 

Remark 5.1. The general case is slightly more involved. It can then happen 
that some eigenvalues of the auxiliary matrix are identical zero or vanish at the 
origin. The relation n = M/2 — S z also ceases to be valid. Particularly notewor- 
thy is the fact that some zeroes might organize in strings (xf, xfq 2 ,..., xfq 2N ~ 2 ) 
signaling degeneracies in the spectrum of the transfer matrix. We refer the reader 
to Ko04c for details and further references. 

Remark 5.2. Note that the identity is the six- vertex analogue of a functional 
relation which has been conjectured in |FM03l [FM04a for the eight-vertex case. 
See the conclusion of Ko04c for further comments and references. The matrix Q + 
is the six-vertex analogue of Baxter's Q-°P er ator i n |Ba72j and Q~ corresponds 
to its counterpart under spin- reversal, compare with the discussion in F M04bj . 
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While this eight- vertex result implies the existence of the six- vertex solutions , 
the algebraic form of Baxter's Q-operator |Ba72) does not allow one to take directly 
the trigonometric limit and obtain an explicit six-vertex Q-operator. 

6. Quantum Wronskian 

In this final section we now wish to make contact with the Bethe ansatz analysis 
of the six-vertex model [L6.7t Su67 by deriving from the results on the spectra 
of the auxiliary matrices (|3.4JI . <|3.13|1 the Bethe ansatz equations. In fact, we will 
actually do more than that: the existence of two linearly independent solutions to 
Baxter's TQ-equation together with (|5.3I) implies a functional equation which leads 
to two sets of Bethe ansatz equations, one above and one below the equator. 

Corollary 6. Let M be odd and N even as before. Renormalizing 

Q±(z) = Q ± (0)- 1 Q ± (z) 

the two solutions of Baxter's TQ -equation then satisfy for any 1 < n < N' 

q nS * Q+(zq 2n )Q-{z) - q- nSZ Q+{z)Q-(zq 2n ) 



(6.1) TW(z) 



qS z _ g S z 



Proof. The identities follow by a simple calculation from l|5.3() and ijOJ. □ 

The case n = 1 of l|6.1|) is special as it only involves the eigenvalues of the 
auxiliary matrices H3.13(l and the quantum determinant H2.12[l . which is explicitly 
known. Hence, this equation, which plays the analogous role of a Wronskian, is 
sufficient to solve the eigenvalue problem of the fusion hierarchy. To stress this 
point we state equation l|6.1|l for n = 1 in a different form. 



e±=e k (xfq 1 ,...,x^ ± q 1 ), k = 0, 1, ^- =F S z 



Corollary 7. For M odd and N even denote by 

M 

T 

the elementary symmetric polynomials in the zeroes of the eigenvalues of ; see 
We use the convention e = 1 and e^ = for k > M/2 =p S z . Then for any 
integer < m < M we have 

/ M X « S*-m+2k_ -S°+m-2k 

( 6 - 2 ) L =E" 



m I f — ' q" — q 

7 fc=o 



Notice that these are M quadratic equations in the M variables {e k } fc=1 

Remark 6.1. Notice that l|6.2() only holds true for an odd number of lattice 
columns M. While we have derived it here for even roots of unity only, it extends 
to the case when q is of odd order as well. However, then the number of solu- 
tions drastically decreases and is in general strictly less than the maximal number 
{m/^-s*) §i ven by the dimension of a fixed spin sector. For example, when q is a 
root of unity of order three there exists only one solution to (|6.2(l ; see Ko04c and 
references therein for details. 

The equations (jfi.2(l do also apply for odd M and q being not a root of unity. 
However, then the representation l|3.1fl defining the auxiliary matrix must be chosen 
infinite-dimensional RW021 IKo04b| . 
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COROLLARY 8. Parametrizing the zeroes x i in as 

(6 ' 3) ~sinf (kf-2n/N) ~ e ^t q ' * 

and setting Bf„ = 1 - 2 cos(^)e i7rk t + e iK ^ +k f" ) we have the identities, 



N 

n± 



M 



(6.4) e" Mfe * [] By = R^" 1 J] B su n ± = Y TS z , 

3=1 i=i 

which are called the six-vertex Bethe ansatz equations above and below the equator. 
In contrast to the equations $6.ty) the Bethe ansatz equations |6'.^| ) are polynomial 
equations of order M in M/2 =p S z variables. 

Proof. Evaluate equation l|6.1|l at n — 1 and z — 1/xfq, 1/xfq^ 1 . □ 



7. Conclusions 

Using representation theory of quantum groups we have derived the Bethe 
ansatz equations of the six- vertex model at even roots of unity and at odd numbers 
of lattice columns M. Using the concept of auxiliary matrices we have shown by 
explicit construction that there are two sets of Bethe ansatz equations and that 
their solutions exist. Moreover, we proved that all the eigenvalues of the fusion 
hierarchy are determined by their solutions. This follows from the fact that for 
even roots of unity and odd M the eigenvalues of the auxiliary matrix (|3.4|1 are 
always non-vanishing polynomials of degree M which must satisfy Baxter's TQ- 
equation. This establishes existence and completeness, albeit we have not shown 
that all eigenvalues of the auxiliary matrix are different, i.e. that there are precisely 
(m/^-s*) different solutions to (j6.4|l respectively Q6.2[l . Our derivation rests on the 
crucial commutation relations i|3.10fl which we proved only for TV = 4, 6 in this arti- 
cle. Nevertheless, in our approach we were able to reduce the question of existence 
and completeness to a well-defined representation theoretic problem, namely the 
existence of an intertwiner for the tensor product of two specific quantum group 
representations. 

We conclude by emphasizing again the concrete practical implications of the 
existence of two linearly independent solutions to the TQ equation. As we derived 
in the text the eigenvalues of Q ± satisfy for the identity l|6.1() which for n = 1 
leads to a quadratic equation in M variables with M being the number of lattice 
columns; see equation (|6.2|l . In contrast, the Bethe ansatz equations l|6.4|l are 
polynomial equations of order M in M/2 variables. For numerical investigations, 
often performed in physical applications, this reduction of the polynomial order of 
the equations is crucial in particular when the system size becomes large M 3> 1. 
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Appendix A. The intertwiners for ir + (z; s) <E> ir + (w; t) when N = 4, 6 

In this section the intertwiners for the tensor product 7r + (z;s) ® ir + (w;t) of 
evaluation representations are stated for roots of unity of order N = 4, 6. The 
construction employs the defining equation of the intertwiner S which is given by 

(A.l) S[n+(z; s) ® ir + {w; t)A(x)] = [{ir+{z; s) ® tt+ (w; t)A op (x)] S, x e 17,(6+) . 

Without loss of generality we can set w = 1 as it follows from the algebraic structure 
of the equations l|A.l|) that the intertwiner S only depends on the ratio of the 
spectral parameters, i.e. S — S(z/w). 

A.l. The case N = 4. For roots of unity of order four the representation 
7r + (z; s) defining the auxiliary matrix is two-dimensional and the equations I A. II) 
are easily solved for the Chevalley-Serre generators. One obtains the result 



S 



( 1 




V o 





l+t-(l+s)Az 
(l + s)(l+z) 
1 + A 
1+z 





Z 





1+A 

1+z 

(l+s)z-(l+t)\ 
(l + s)(l+z) 





\ 




A / 



Here A is a free parameter and the matrix is computed with respect to the basis 
{|0,0} , |0, 1} , 1 1,0) , |1, 1)} in the tensor product ir + (z;s) <g n + (l;t). The arbitrary 
overall normalization factor of S is fixed by setting the matrix element of the highest 
weight vector equal to one. 



A. 2. The case N = 6. Now the representations tt + (z; s), tt + (w — l;t) are 
each three-dimensional. As S commutes with the Cartan generators we immediately 
deduce that S must be block-diagonal w.r.t. the following decomposition of the 
tensor product space, tt + (z;s) (8 ir + (l;t) — Vi © V2 © V3 where the respective 
subspaces are spanned by the following basis vectors 

Vi = span{|0,0),|l,2),|2,l}}, 
V 2 = span{|0,l),|l,0},|2,2)}, 
V 3 = span{|0,2),|l,l},|2,0)} . 

The calculation is cumbersome but straightforward and one finds the following 
solution up to a common normalization factor, 



S\ 



s\ 



s\ 



v 3 



1 



(1- 





(l-z)(zq- 1 +t) 
(sz+q)(sz+q- 1 ) 

z{zq+t)(t+q) 
(sz+q)(sz+q^ 1 ) 

z)(szq- 1 + l 



(sz+q)(sz+q- 
z(sz+q)(t+q" 
(sz+q)(sz+q- 

V 

f (l-z)(l+zq) 
(sz+q)(sz+q 1 ) 

(t+q)(l-z)z 
(sz+q)(sz+q- 1 ) 
{t+q)(t+q- 1 )z i 
(sz+q)(sz+q- 1 ) 



\ 

(q+s)(t+zq~ 1 ) 

(sz+q)(sz+q- 1 ) 

(sz-t)(z+tq) - 
(sz+q)(sz+q- 1 ) / 

(sq+^jl+szq- 1 ) 
(sz+q)(sz+q- 1 ) 
{sz — t)(sz+q) 



(sz+q)(az+q~ 





(zq+t)(zq M-t) 
(sz+q)(sz+q- 1 ) 



q(w-l)(sq+l) 



{sz+q)(sz+q 
w{l+s)(l+t)-tq- L - 
(sz+q)(sz+q~ 1 ) 
(t+q- 1 )(sz-t)z 
(sz+q)(sz+q- 1 ) 



w~ sq 



(s+q^s+q- 1 ) 
(sz+q)(sz+q~ 1 ) 

q(s+q)(sz-t) 
(sz+q)(sz+q~ 1 J 
(sz — t) (sz+tq 



(sz+<j)(sz+ij- 
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